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Abstract 

A method to calculate particle fluxes applicable in most of the spectroscopy techniques is de- 
scribed. Flux intensities of backscattered or absorbed electrons and emitted photons are calculated 
using a method of convergence to solve the Invariant Embedding equations that are used to describe 
the particle trajectories inside a solid sample. Our results are found to be helpful to carry out a 
procedure for quantitative characterization using instrument such as Electron Probe Microanalysis 
or other probes. Examples of application to calculate the composition of ternary alloys are given 
and are compared with the same calculations using another procedure. 

PACS numbers: 73.20.Mf; 72.20.-i 
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I. INTRODUCTION 



In most of the spectroscopy techniques, the signals produced by the different probes must 
be recorded and interpreted considering the different kinds of interactions between particles 
during its trajectory through an analyzed sample. This procedure involves hard calculations 
of fluxes of particles like electrons and/or photons, neutrons, protons etc. These fluxes 
must be obtained solving a transport problem where some particles produce other particles: 
as an example, electrons transport could produce electrons, photons or both of them as a 
consequence of their interactions with matter. The transport problem, in the spectroscopy 
area, is usually solved using Monte Carlo Calculation |lH4| or numerical solutions of the 
the time-independent Boltzmann transport equation. Although considerably advances were 
made in the subject, most of the methods used in the practice have an empirical or semi 
empirical origin |5l-ll0]. 



The quantitative determination of the unknown chemical composition of a material is an 
essential goal for its characterization. Among the spectroscopy techniques, the most used 
by researchers in materials science for compositional determinations is the Electron Probe 
Microanalysis (EPMA). EPMA is based on the detection of the characteristic X-rays emitted 
from a solid sample when a beam of electrons impinges on its surface with energies that 
typically vary between 2 and 50 keV. This instrument combines the capabilities of both the 
scanning electron microscope and an X-ray fluorescence spectrometer. The electrons of the 
beam are scattered as they penetrate within the sample producing backscattered electrons, 
secondary electrons, absorbed electrons and characteristic and continuum X-rays spectra. 
Quantitative microanalysis in EPMA is performed by comparing the characteristic X-ray 
intensity from one element in the sample with that from a reference standard containing a 
known amount of the same element, by means of the so called the K ratio, 

K = T (1) 
the characteristic X-rays emitted from the analyzed 



Where I and I p are the intensities o 



sample and that from a standard. [2j, lll|. 



12 1 . The K ratio could be transformed into the 



mass concentration of the element of interest taking into account some correction factors 
due to processes undergone by the characteristic X-ray and by the electrons of the beam 
inside the sample. These processes involve scattering dynamics related to the sizes of the 
different atomic nuclei present in the sample and the standard, absorption, fluorescence. 



The respective correction factors arising from each of these processes are called Z,A and 
F in the literature [2(. Then, the atomic fractions of the elements are obtained using the 
following expression: 

K = f = ZAF^ (2) 

Where C and C p are the atomic fractions of the element of interest in the sample and in the 
standard. 

The most common procedures to estimate the corrections are the so called Matrix Cor- 
rection which are based on the calculation of the X-ray production as a function of the 



sam 
Q,[l3| 



e depth, tp(pz) to estimate both the absorption and the atomic number corrections 



14|. Most of the analytical expressions for the <f{pz) functions used in commercial 



software packages are obtained from ad hoc models. Also the expressions for the ip(pz) 



15|. On 



function and for the backscattering electron coefficients are consider as unrelated 
the other side, in the literature there are theoretical models based on Boltzmann's trans- 
port elation to describe the phenomenon and to intent the signals in EPMA Q Q. 
Although, this procedure has the disadvantage that there are solutions only for the one di- 
mension problem, and little progress has been made on the 2-D or 3-D cases. Therefore this 
has a reduced practical application. Monte Carlo calculations is still the most worthwhile 
method of approaching theoretical problems in EPMA j^, 12, Q, 19]. However, depending 



on the accuracy or resolution required, running times can be computationally expensive. 

In the last two decades a few authors demonstrated the efficiency of the Invariant 
bedding (IE) method to describe the scattering of electrons inside a solid sample 
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27|. 



This method is useful to estimate the flux of particles through the interfaces in a material 
medium as a function of experimental parameters and the detected signals. It is, the IE 
method avoids calculations of fluxes inside the analyzed samples, thus gaining calculation 
advantages, and yields a set of equations for the detected fluxes as a function of the size 
of the analyzed sample. This set can be solved as an initial- value problem. Particularly, 
;he backscattering problem was successfully solved by different authors using IE method 



22j, |28(. However, a procedure based on IE, useful for the experimentalist to perform chem- 
ical characterization, is still lacking. 

In this article we focus on the calculation of the chemical compositions of a materials 
by using the experimental K ratios and theoretical X-ray intensities obtained with IE. Al- 
though the IE method is a powerful tool in order to simplify the transport equations, in 



3-D transport problems the addition of degrees of freedom of the particles generates more 
complex mathematic. The Method of Convergence (MC) proposed here allows us to deduce 
closed approximated solutions for the 3-D IE equations for the detected X-rays intensities. 
These intensities can be calculated with no need of previous knowledge of the (p(pz) function 
because the absorption corrections are already included in the equations. Here we propose 
two solutions for the IE equations with different levels of approximation. One called the 
simple expression and the other, more rigorous, called the complete expression. Compari- 
son with other model proves that both solutions are a suitable tool for the calculation of 
intensities and concentrations. 

II. THEORETICAL MODEL 

The states ladder model and MC were are described in previous reports[25(]. Basically, 
after impinging with energy EO, electrons can suffer two kinds of interactions: elastic or 
inelastic. Elastic collisions only cause deviations while inelastic ones only produce energy 
losses. Based on the assumption that it is allowed only one kind of inelastic collision with a 
constant energy loss of value AE, a discrete set of possible values is generated: El, E2, E3 and 
E4, for a five level system. For energies E < E4, electrons are deactivated. Characteristic 
radiation can be generated exclusively by the inelastic collisions. Contributions 101, 102, 
103, 104 to the total intensity are made by electrons backscattered with El, E2, E3 and 
E4, while contribution 14 corresponds to the absorbed electrons. All these intensities can 
be calculated with the IE method. This method starts with the deduction of the integro- 
differential equation for each contribution to the intensity, whose derivation variable is the 
thickness r of the sample. In the case of / Q o,^n the intensity contribution for the electrons 
that impinge with EO and «o and are backscattered with El and fli, the following equation 
is obtained: 
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, flO^Wo^K , 2( g0 + (T ) 



I 



«0 7e e l 
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e l 7e e 




(3) 



/ 7 J 8 71 

Were the simplified notation dj = cos(aj) and (3j = cos (j3j) is used for the incident and 
emerging director cosines of the polar angles of the electrons, measured with respect to the 
normal of the the sample surface. The index j = 0, 1, .. labels the energy states values, <Tj 
and Sj are the elastic and inelastic cross sections, respectively; ip, e and 5 are the scattering 
angles due to single interactions inside the sample, and gj is the electrons efficiency of X-ray 
production per collision. 

Using the MC it is possible to obtain two approximate solutions for the last equation: the 
complete form, i.e., the formal solution, and the simple form, used initially as an auxiliary 
expression to solve the equations of the lower levels. These approximate solutions tend to 
an exact solution for r — > oo.The simple form is: 

/ \ god /3is (i + A)ri 
l d a (oo) = — (4) 

(A + l)(si«o + s $i) + /jidoPi Mo 



Where J = A + 2 - aTi, r\ = A = 2 - 2 V / ATT and A = a /s . 

In order to study the applicability range of the solutions, chemical composition calcu- 
lations made with the simple and the complete expressions are presented. The factor of 
generation gj it calculated starting from Bethe's expression for Lk lines, and an energy level 
Ej, for the ionized element in the compound: 



. p k N A E s 2E S 



Where k is the partial density , A^is the atomic number and m is the mass of the electron. 
Parameters such as the partial densities, ionization energies or attenuation coefficients are 
obtained from tables. The cross sections are calculated from classic expressions as functions 



of the partial densities. In this model, it is considered that the beam electrons interact with 
atoms in the compound, whose cross sections is the pondered average of the cross sections 
of each element. The procedure is as follows: 1) The summation of the measured K ratios 
is normalized. 2) This normalized values are considered as the initial atomic concentrations. 

3) With these new initial atomic fractions, the K ratios are recalculated with the simple 
and the complete MC expressions, and the results are compared with the measured values. 

4) In successive calculations, the atomic fractions of both expressions are adjusted in order 
to fit the theoretical and experimental K ratios. The final values for the concentrations are 
those that yield the best fit. The calculation was performed for the line L a of Zr and/T a of 
Fe. 



III. CHARACTERISTIC INTENSITIES. SIMPLE EXPRESSIONS. 

We proceed to calculate 101, 102, 103, 104 and 14 in order to obtain / = ^2 { I0i 
The simple expressions for the Intensities are: 

701 (oo) = t 4floa 5o/3i(l + A)Tid/3i 

Jo {(A + l)(sid? + s /3i) + /idoMAJo 

J02(oo) = I' %2a g (l + X)T 1 W $ 2 d$ 2 

Jo [(s 2 d + s /3 2 )(l + A) + /idofyXJo 



(8) 



703(oo) = 

1 16g 3 a s (l + X)Tj 3 d$ 3 {W Ro + b c Z (\ + 2)} 
(s 3 d + s f3 3 )(l + A) + fid f3 3 \J^ 

/04(oo) = f 1 32 W oT 1 (A + l)/3 4 
Jo {AJ^[s 4 cfo + So/34]} 
x{b c J Z (d V + A + 2) + + J ) + 

aMiWoRo + coZ [b (\ + 2) + d V + W ]}}d^ (9) 

where 



_ a (6o-2)T 1 + 2A + 4 

'° A + 2 + 2\/ATT - a T! 
^ = (2 + A)(6 + l)-a Ti 



(10) 
(11) 
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IV. ASSESSMENT OF THE RESULTS: ZR-BASED ALLOYS 



In order to compare the calculated K ratios and concentrations between MC and PAP 



29- 



311 ] results, we used experimental data obtained with an analytical microprobe of electrons 



CAMECA SX50 with 20 kV of accelerating potential from samples of ternary metallic alloys 
of the systems Zr-Sn-Fe and Zr-Sn-Nb. For the calibration, high purity standards of each 
of the analyzed elements were used. The proposed method was applied to a set of fourteen 
measurements on these alloys. The nominal atomic composition of the Zr-Sn-Fe sample 
is: Zr 62.5 - Sn 25.0 - Fe 12.5. It was thermally treated to 800C for 2000 hours. Three 
phases are detected in this sample: two binary intermetallic compounds (ZrFe2 and Zr 5 Sn 3 ) 
and a ternary compound Zr 6 Sn 2 Fe. At first, the atomic composition was estimated using 
the PAP procedure. On the other hand, the theoretical results of section II and III were 
used to estimate the theoretical K ratios and the atomic composition with the complete 
and the simple expressions for the characteristic intensities. The results are shown in tables 
I and II. They show comparisons between PAP concentrations and results obtained with 
the simple and the complete 3D-IE expressions, for the Zr 5 Sn 3 and Zr 6 Sn 2 Fe phases. For 
each: element of interest, measured K ratio (k), concentrations from the PAP method 
(at (PAP)), concentrations from 3D-IE simple equation (at (3Ds)), concentrations from 3D- 
IE complete equation (at(3Dc)), K ratios from 3D-IE simple equation (ks), K ratios from 
3D-IE complete equation (kc). Each table contains different sets of values, corresponding to 
different measuring points in the same phase. The first rows of each set of values indicate 
the PAP, 3Ds and 3Dc estimations of the oxygen concentration, as this elements is always 
present at the surface of these alloys. 

On the other hand, the nominal atomic compositions of the Zr-Sn-Nb samples are: Zr 
70.0 - Sn 5.0 - Nb 25.0 for the first one, and Zr 65.0 - Sn 10.0 - Nb 25.0 for the other one. 
They were thermally treated to 950C for 2900 hours. Table III shows comparisons of the 
theoretical results in this article and experimental data, in the same way as Tables I and II. 



7 



TABLE I: Theoretical calculations (IE) for the K ratios of the elements (E) in the ZrsSn3 phase of 
the Zr-Sn-Fe alloy. The 3D-IE results obtained with the complete (c) and the simple (s) expression 
are compared with PAP calculation. 







3.16% O 


2.6% O 


15.1% O 








K 


%at(PAP) 


%at(3Ds) 


% at (3Dc) 


K s (t) 


K c (t) 


Zr 


0.4498 


61.3297 


53.1501 


60.9755 


0.4460 


0.4791 


Sn 


0.3303 


38.2826 


46.4461 


38.4317 


0.3268 


0.3418 


Fe 


0.0020 


0.378 


0.4038 


0.5928 


0.0020 


0.0020 






6.28% O 


1.1% O 


13.1% O 






Zr 


0.4527 


61.7237 


53.93 


61.4206 


0.4573 


0.4883 


Sn 


0.3336 


37.9375 


45.6761 


38.0250 


0.3243 


0.3486 


Fe 


0.0018 


0.3387 


0.3557 


0.5544 


0.0018 


0.0019 



V. DISCUSSION 



Calculations performed with the simple and the complete expressions obtained with MC 
are compared with results of the PAP method. The fitting is acceptable keeping in mind 
the conditions of calculation and the fact that the intensity expressions were deduced from 
approximate expressions for the transport coefficients of the electrons. Thereby, the inten- 
sities are a second stage of approximation. In the Zr-Sn-Fe results, the presence of oxygen 
changes the composition of the phases in a significant way. Oxygen diminishes the K ratios 
of the other elements, in comparison with non-oxidized samples. In relation to the phase 
Zr 5 Sn 3 , it is necessary to remark that the results yielded by the complete MC expressions 
fit the PAP concentrations better than the simple ones. However, the complete expressions 
overestimate the oxygen concentrations at the surface. In relation to the MC calculated 
K ratios, the discrepancy is less than 5%, with respect to the PAP values. In general, for 
this compound, the simple MC expressions underestimate the Zr and overestimate the Sn, 
in comparison with PAP and complete MC expressions. It is necessary to emphasize that 
the equations are valid if the ratio ip = fi/s can be considered as negligible with respect 
to unity. In the calculations,^/ values are: 0.05% for Fe, 0.18% for Zr and 0.14% for Sn. It 
means that the MC expressions are at the limit of it applicability for the present selection 
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TABLE II: Theoretical calculations (3D-IE) for the K ratios of the elements (E) and its atomic 
fractions in the Z^S^Fe phase of the Zr-Sn-Fe alloy. The 3D-IE results are compared with PAP 
calculation. 







1.14% O 


7.0% O 


7.0% O 








K 


%at(PAP) 


%at(3Ds) 


% at (3Dc) 


K s (t) 


K c (t) 


Zr 


0.5716 


67.36 


61.2063 


65.294 


0.5561 


0.5982 


Sn 


0.2120 


22.7 


27.7656 


21.6365 


0.2090 


0.2205 


Fe 


0.0604 


9.925 


11.028 


13.071 


0.0605 


0.0625 







1.47% O 


5.5% O 


7.0% O 






Zr 


0.5763 


67.43 


62.1693 


65.3326 


0.5692 


0.5990 


Sn 


0.2115 


22.5 


27.2487 


21.5054 


0.2065 


0.2193 


Fe 


0.0616 


10.05 


10.5820 


13.1720 


0.0585 


0.0631 







1.18%0 


7.4%0 


15.4% 






Zr 


0.5102 


65.0721 


58.1445 


63.4064 


0.4985 


0.5228 


Sn 


0.2680 


30.1344 


36.6928 


29.8788 


0.2621 


0.2797 


Fe 


0.0272 


34.8033 


5.16274 


6.7148 


0.0268 


0.0285 







2.13%0 


4.0%O 


7.5% 






Zr 


0.5625 


66.5781 


60.6771 


64.5946 


0.5561 


0.5870 


Sn 


0.2230 


23.7969 


29.1667 


22.7027 


0.2217 


0.2304 


Fe 


0.0585 


9.6352 


10.1562 


12.7027 


0.0562 


0.0603 







0.21%O 


2.3% O 


9.0% 






Zr 


0.5629 


66.4305 


60.0338 


64.4540 


0.5783 


0.5784 


Sn 


0.2292 


24.3489 


29.9064 


23.3171 


0.2340 


0.2341 


Fe 


0.0562 


9.2206 


10.0596 


12.2288 


0.05729 


0.0573 
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TABLE III: Theoretical calculations for K a ratios of elements (E) in the ZrSnNb alloy. . 





K 


%at(PAP) 


%at(3Ds) 


% at (3Dc) 


K s (t) 


K c (t) 


Zr 


0.6071 


64.50 


58.92 


59.02 


0.6129 


0.5946 


Sn 


0.0220 


2.59 


1.74 


2.44 


0.0217 


0.0228 


Nb 


0.3253 


32.91 


39.34 


38.54 


0.3275 


0.3155 



Zr 0.616 66.03 63.62 61.32 0.6337 0.6034 

Sn 0.0678 7.81 5.72 6.74 0.0696 0.0671 

Nb 0.2552 26.17 31.66 31.94 0.2556 0.2529 



Zr 


0.5957 


64.88 


59.02 


59.02 


0.6126 


0.5934 


Sn 


0.0271 


3.26 


2.16 


2.84 


0.0269 


0.0269 


Nb 


0.3076 


31.86 


38.82 


38.14 


0.3224 


0.3115 



Sample 2 



Zr 


0.7291 


78.39 


74.77 


74.77 


0.7464 


0.7309 


Sn 


0.0380 


4.46 


3.21 


4.09 


0.0325 


0.0382 


Nb 


0.1677 


17.14 


22.02 


21.14 


0.1746 


0.1629 



Zr 0.6931 70.82 67.01 67.01 0.6831 0.6661 

Sn 0.0294 3.51 2.99 2.99 0.0292 0.0278 

Nb 0.2489 25.67 30.00 30.0 0.2487 0.2390 

Zr 0.7071 73.64 68.78 68.10 0.6993 0.6762 

Sn 0.0269 3.10 2.20 2.90 0.0268 0.0269 

Nb 0.2344 23.26 29.02 29.0 0.2351 0.2304 



Zr 


0.6849 


72.92 


68.10 


68.10 


0.6932 


0.6762 


Sn 


0.0277 


3.25 


2.25 


2.90 


0.0274 


0.0269 


Nb 


0.2352 


23.255 


29.65 


29.0 


0.2405 


0.2304 
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of parameters. In relation to the oxygen, the simple MC expressions overestimate its con- 
centration in comparison with PAP. The results for Zr 6 Sn 2 Fe show the same behavior that 
those of the previous phase in the sense that there is a better fit between concentrations 
calculated with the complete MC expressions and PAP in comparison with the simple ones. 
These expressions yield results for the K ratios that differ from the experimental values in 
an amount that ranges from 2 to 5 %. In general, the theoretical K ratios fit to experi- 
mental values within less than 3 %, with an overestimation in the results of the simple MC 
expressions. The ip ratios for this phase are: 0.05 for Fe, 0.17 for Zr and 0.14 for Sn, so the 
MC expressions are at its limit of application. With regard to oxygen, PAP concentrations 
are about 1% while MC values range from 7 to 15%. In this phase the discrepancies in the 
concentrations of Zr, Sn and Fe between PAP and MC range from 2 to 5%. In the case 
of Zr-Sn-Nb, the fit between K ratios is better, but the discrepancy between PAP and MC 
concentrations is higher. Both MC expressions underestimate the Zr and overestimate the 
Nb in comparison with PAP results. There is no oxygen at the surface of this compound, 
and the \l/ ratios are: 0.12% for Nb, 0.16% for Zr and 0.17% for Sn. It is important to remark 
that the inelastic-elastic correlation determines the ionization rate and the <p(pz) function. 
Thereby, the intensities depend on such correlation, which is related to the A factor, this is 
to say, to the ratio between the inelastic and elastic cross sections. Our calculations were 
performed basically to verify the internal coherence of the 3D model and the MC expressions. 
This means that parameters such as cross-sections, attenuation coefficients, ionization ener- 
gies, etc are fixed values taken from tables, and the only adjustable parameter is the Emin 
value, i.e., the minimum energy lost per inelastic collision, whose magnitude was selected to 
fit the backscattering and absorption coefficients of the electrons. Bearing in mind this fact, 
the presented results for concentrations and K ratios are quite satisfactory and indicate that 
the 3D model and the MC expressions are appropriate for the theoretical representation of 
the interaction phenomena in EPMA, and thereby, for the characterization of compounds 
by means of this or similar techniques. 
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Appendix A: Determination of the function I ao ^(oo). Method of Convergence. 



Considering in the integral-differential equation for /^^(oo), Eq.Q the six first terms 
equals to a constant II: 

g R00 a ps 



II 



0i 



"0 Jtp «0 



c,Jn Co 



Eq.fl3|) is now written as 



Bl/3 

+ <Tl) + 

Equation [Al] contains four integrals: 



(Al) 



— =n-/ a0jj8l (r)(— ; + — ^ + //) (A2 

dr a fi x 

The solution for r — >• oo is: 

4 A (oo) = - - (A3) 

«o(si + cti) + A*i(s + a ) + yuaoPi 

Let us rename the term II x cf =B1, were Bl=Bl(a , 



I aoA (oo) = ^ (A4) 

«o(Sl + CTi) + /3i(Sq + (T ) + /i«0/3l 



Q'o 

(A6) 




ic A Co 

These integrals are solved by replacing / QOj/3l (oo) by the expression [A4] obtained with 
the MC. The functions of Elastic Backscattering R00 are replaced by its normalized simple 
expression, deduced with the MC [reference]. Replacing at [Al], solving, multiplying for «o 
and replacing — > ty — /V s o an d ^1 = /V s i we obtain. 
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2g a soT 1 XsqBIPx 
dV = h 



A 2si[A + l + *i/9i] 

x ln( f 1 + 1 + ^° ) + aoSi^lTiaoA 



^is ' A + r s [A + l + ^o«o] 
x ln(- h 1 + -) ln(l + — ) + 



"osi A + l ft si[A+ l + *i4i] 

xln(^- + l + -^-)ln(l + -i) 

Asifildo . t s *i . 

r /n(- h 1 + -) 



2s [A + 1 + ^ «o] «oSi A + l 

Substituting this results in [A4] we obtain the following expression for 101: 

k 



/01 Q/3 (oo) c 



ao(si + + &(so + °o) + A*aoA 



25f aos Ti s A/3i Si 
I ; h Aso^ofi'ol — ln(- — 

+ 1 + TTT) + I rx I i I , T , i i ln (-T— + 1 + 



A + l 7 s [A + l + ^ «o] «oSi A + r 

wi , ^ \ i aoSoTiao/^i , si ^ v 
x ln(l + — ) H —ln(- h 1 + -) 

/V Sl [A + l + *Mi] /3is A + l 

x ln(l + -) + „ Agl "° + 1 + A)} 

a 2s [A + 1 + ^ «o] a si A + l /J 



(A9) 



Where the super index "c" it refer to "complete". 

Let us now substitute Bl = K Q g Q \s at [A4] and a simple expression of 101 is written as: 

rm / \s K g \s Q (3 1 

i01 Q/9 (oo) = — 

«o(si + + A( s o + cr ) + fia (3i 
Where the super index s it refer to simple. In order to determine kq we use the equality: 

/ / I01 a ^oo) s da dk = [ [ J01 a/3 (oo) c AM/3i 
Jo Jo Jo Jo 

Substituting [A9] and [A10] in [All], considering s = «i and /i « Sj, 
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' ■" 1 k g \ da dfii f 1 f 1 g /3i 



J (A + l)(a + /3i) Jo Jo (A + l)(& + /3i) 

2a T! r a ^itto/3i, n 1 x, / 1 , V \ 

{ _ + « A{- rn -Zn(l + £ )Zn(^ + l) 



A + 1 ^ d ' 

M ° -H- + l) + 



2(A + 1) v ft 
Solving the integrals k is written as: 

4a (A + 1)7\ 



2(A + 1) v d 
ln{- — h 1)}}c?q;o 



Kq 



A 2 Jo 



(A10) 



(All) 



Substituting [A13] in [A9] and [A10], the simple and complete expressions for 101 are 
obtained respectively. 
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